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ABSTRACT 



Aims. We study the chemical complexity towards the central parts of the starburst galaxy M 82, and investigate the role of certain 
molecules as tracers of the physical processes in the galaxy circumnuclear region. 

Methods. We carried out a spectral line survey with the IRAM-30m telescope towards the northeastern molecular lobe of M 82. It cov- 
ers the frequency range between 129.8 GHz and 175.0 GHz in the 2 mm atmospheric window, and between 241.0 GHz and 260.0 GHz 
in the 1.3 mm atmospheric window. 

Results. Sixty-nine spectral features corresponding to 18 different molecular species are identified. In addition, three hydrogen recom- 
bination lines are detected. The species NO, HiS, H2CS, NH2CN, and CH3CN are detected for the first time in this galaxy. Assuming 
local thermodynamic equilibrium, we determine the column densities of all the detected molecules. We also calculated upper limits 
to the column densities of fourteen other important, but undetected, molecules, such as SiO, HNCO, or OCS. We compare the chem- 
ical composition of the two starburst galaxies M82 and NGC253. This comparison enables us to establish the chemical differences 
between the products of the strong photon-dominated regions (PDRs) driving the heating in M 82, and the large-scale shocks that 
influence the properties of the molecular clouds in the nucleus of NGC 253. 

Conclusions. Overall, both sources have different chemical compositions. Some key molecules highlight the different physical pro- 
cesses dominating both central regions. Examples include CH3CCH, C-C3H2, or CO + , the abundances of which are clearly higher in 
M 82 than in NGC 253, pointing at photodissociating regions. On the other hand, species such as CH 2 NH, NS, SiO, and HOCO + have 
abundances of up to one order of magnitude higher in NGC 253 than in M 82. 
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1. Introduction 
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Despite its small size, M82 is the bright est IRAS galax 
(Sim = 135 1.02 Jy; L m = 3 x 10 10 L G ; iRice et alJ TT98~ 
Telesco, 1988) . Its cent ral starburst region, of ~ 500pc size 



dO'Connell & M angano 1978), has abolometric luminosity 100 
times higher than a region of equiva l ent size in our Galax y 
dRieke et aU 119801: iLester et alJ 119901: lO'Connell et ail Il995h . 
Almost all its luminosity arises from the central kiloparsec 
(lLo etalJll987h . which seems to be filled with a large number 
of compact star clusters whe r e high-mass stars dominate the 
luminosity dRieke et all Il980t iKnapp et alJ fl980l: iFuente et all 
2008; Westmoauet teet al.ll2009t) . The galaxy M82 is consid- 
ered a prototype starburst galaxy because o f its high star-forming 
rate (~ 9 Mm yr -1 , IStrickland et alj|2004l) . which was probably 
caused by a recent interaction w ith its companion M 8 1 , at a 
projected angular distance of 36' ( Cottrelll[T977b iBrouillet et alJ 
Il99ll ISofudl 1998b ISun et alJl2005b . The M 82 burst of star for- 
mation is thought to be fuelled by large am ounts of molecular 
gas (M Hz ~ 2 x 10 8 Mm. iNavlor et aljfeoiOh still present in the 
central 1 kpc region. This region shows two emission peaks op- 
posite the dynamical centre, the north-east (NE) and south-west 



(SW) molecular l obes dMaoetalJ 120001: iGarcfa-Burillo et al.l 
I2002t lMartfnetalJ l2006a). which are likely part of a nuclear 
ring of warm ga s and dust heated by stars dNakai et alj|1987h 
iDietz et alJ[l989T) . 

Apart from its enhanced activity, the proximity of M82 
(D = 3.6 Mpc, iFreedman et alJ fl994l) ensures that it is one 
of the extragalactic sources studied in most detail at all 
wavelengths, from radio to X-rays. M82 shows a bipo- 
lar asymmetric out flow of cool (atomic and molecular), 
warm and hot gas dLvnds & Sandagel 1 1963b iBreg man et al. 
U995blShopbell & Bland-Hawthornl 19981) . which extends above 
and below th e galactic plane to large scales, even reach- 
ing the halo dSeaauist & Odegardlll99lb lMcKeithetal.ll 1995b 
ISeaquist & ClarkfeOOlblVeilleux et al.l l2009). and creating chim- 
neys, supershells, and bubbles of atomic and molecular material 
dGarcfa-Buriflo et al.ll2001blWills et alJ20"02h . These superwinds 
are generated by the nuclear starburst, which, accompanied by 
a high supernova rate, creates large photon-dominated regions 
(PDRs) in the central parts of the galaxy. At the same time, this 
galaxy is a strong emi tter of X-rays (e.g.lRead & Stev ens 2002) 
and cosmic -rays (e.g. IVERITAS Collaboration et alJt2009l) . al- 
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though they do not seem to domin ate the chemistry of the M 82 
nucleus as much as the UV field dFuente et al.ll2008h . 

The UV field in the nuclear region of M82 character- 
izes the chemistry of its molecular gas. Many species show 
lower abundances in M 82 than in other nearby starburst galax- 
ies (such as NGC253, NGC4945, or NGC6946). In par- 
ticular, the more complex ones (with > 3 atoms) such as 
CH3OH, NH3, or HNCO, seem to be underabun dant, since 
they a re easil y dissociated by the strong UV fields dRieke et al.1 
1980. 119881: iNguven-O-Rieu et alJ mm. iTakano et all 12002: 
Mauersberger e t al. 2003; M artin et al.ll2009al) . In particular, the 
HNCO abundance shows large variatio ns of nearly two ord ers of 
magnitude among starburst galaxies (Martin et al. 2009a). This 
suggests that the observed HNCO abundances are related to the 
evolutionary stage of their nuclear starbursts. 

The chemistry of the M 82 nucleus seems to be the result 
of an old starburst mainly dominated by the photod issociating 
radiat ion from the already-formed massive stars dFuente et al.l 
120081) . Thus, most of the dense molecular material for fu- 
ture sta r-forming events ha s already been consumed in this 
galaxy ( Marti net al.l 1200 6a). Other galaxies, notably the ones 
cited above, are interpreted as being at an earlier stage of a 
starburst, where PDRs, are significant (Martin et al. 2009b), 
but do not yet play a major role. They are instead domi- 
nated by large-scale shocks among molecular cloud complexes, 
caused by cloud-cloud collisions and ma s s loss from mas- 
sive stars dGarcia-BuriUo et all l2000L 12001 IWang etaD l2004t 
lAladroetal.ll2011al) . 

To confirm that the evolutionary states of starburst galaxies 
can be characterized by their chemical compositions, we per- 
formed 2 mm and 1.3 mm spectral scans towards the NE molec- 
ular lobe of M 82. In Sect. [2] we present the observations and 
the data reduction. The data analysis and the estimated uncer- 
tainties are explained in Sect. [3] In Sect. [4] we present the de- 
tected molecules as well as some important undetected species. 
Section [5] shows a detailed comparison of the chemistry that we 
found in M 82 and NGC 253, using a surv ey of similar f r equenc y 
performed towards the second galaxy by iMartin et al.l d2006b). 
Finally, in Sects. [6] and [7] we draw our conclusions and sum- 
mary. 

2. Observations and data reduction 

Observations were made with the IRAM 30-m telescop43, lo- 
cated at Pico Veleta, Spain, during several periods between 2006 
and 2009. The frequency survey was carried out towards the NE 
molecular lobe of the starburst galaxy M 82, which has an offset 
position of (Aa, AS) = (+13.0", +7.5") with respect to the dy- 
namical centre (a 2 ooo = 9 A :55 m :51 J .9Q, 5 2 ooo = + 69°:40':47".10; 
iJov et al.l ll987) (Fig. Q]). At this molecular cloud complex, the 
dissociating radiation is clai med to be stronger than i n the SW 
lobe and the centre regions (iGarcfa-Burillo et aDl2002l) . We cov- 
ered the frequency range from 129.8 GHz to 175.0 GHz at 2 mm 
wavelengths, and from 241.0GHz to 260.0GHz at 1.3mm. We 
simultaneously used the (now decommissioned) C and D re- 
ceivers in single sideband mode, which provided an instanta- 
neous bandwidth of 2 x 1 GHz at both 2 and 1.3 mm. As a 
backend, we used the 256 x 4 MHz filterbanks. The observations 
were made by wobbling the secondary mirror with a switch- 
ing frequency of 0.5 Hz and a beam throw of 220" in azimuth. 
The data were calibrated using the standard dual load method. 

1 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



The rejection of the image sideband was measured for each fre- 
quency tuned, and then varied from 9dB to 24 dB. The pointing 
was checked every one or two hours by observing planets and 
nearby bright continuum sources. We estimate the pointing accu- 
racy to be ~ 3". The focus was always checked at the beginning 
of each run and during sunsets. The beam sizes ranged from 19" 
to 14" in the 2 mm band and from 10" to 9" in the 1.3 mm band. 
The spatial resolution ranged between 158 kpc and 333 kpc. 

The observed spectra were converted from antenna temper- 
atures (T^) to main beam temperatures (7mb) using the rela- 
tion Tmb = (F e ff/B e s) T* A , where F e ff is the forward efficiency 
of the telescope (which was 0.93 and 0.88 for the 2 mm and 
1 .3 mm bands, respectively) and B e ff is the main beam efficiency, 
ranging from 0.72 to 0.64 for the 2 mm band, and from 0.50 
to 0.43 for the 1.3 mm band. Figure [2] shows the reduced data. 
Baselines of order or 1 were subtracted in most cases. We 
achieved an rms of ~ (2 - 5) mK across most of the 2 mm band, 
and of ~ (3 - 8)mK for the 1.3 mm band (for the above men- 
tioned 4 MHz channels corresponding to ~ (7 - 9)kms~ 1 and 
~ (4 - 5)kms _1 , respectively). 

3. Data analysis and uncertainties estimates 

Assuming local thermodynamic equilibrium (LTE) and optically 
thin emission for all the detected species, we plotted Boltzmann 
diagrams in order to estimate the rotational temperatures (r rot ) 
and the molecular colum n densities (Af mo i). The deta ils of the 
method are explained in iGoldsmith & Lang er (1999). The re- 
sulting rotational diagrams are shown in Fig. lA.2l and the derived 
physical parameters are listed in TableQ] 

For the LTE analysis, the data were also corrected for beam 
dilution by Tg = [((% + 0?) / 0^] Tmb, where T% is the source- 
averaged brightness temperature, 9 S is the source size, and 6\, is 
the beam size in arc-seconds. 

3. 1 . Opacity of the detected molecules 

The rotational diagram method assumes optically thin lines. To 
conclusively prove that extragalactic lines are actually optically 
thin is, however, not an easy task because the ratios of hyperfine 
structure lines can be used only in those exceptional cases where 
the separation between these line s is larger than the Doppler 
width. This is the case for CN (e.g. lHenkel et alJll998l) and also 
for NO, but not for most other molecules with hyperfine struc- 
ture. A second possibility is to compare line intensity ratios with 
known isotopic ratios: if there is no isotopic fractionation be- 
cause of chemical effects (which are basically expected to occur 
at lower temperatures) or selected photodissociation, a line in- 
tensity ratio of the same transition of two isotopologues that is 
close to the isotope ratio can be taken as an indication that both 
lines are optically thin. The problem for extragalactic studies is: 

a) that isotopic lines (in particular for molecules that contain two 
substitutions of rare isotopes) may be too weak to be detected; 

b) the isotopic ratios are not a priori kn own and may differ from 
those expected in our own Galaxy (e.g. lMartfn et al.l l2010): and 

c) the sources are not resolved, hence the opacity cannot be de- 
termined from the brightness temperature and the excitation tem- 
perature. The lines of 12 C 16 are the only millimeter lines for 
which we can be sure that these are optically thick in many ex- 
tragalactic sources. The fact that all the other mm-lines in extra- 
galactic sources tend to be more than a factor of ten weaker than 
CO may be an indication, but no proof, of optical thinness, since 
the low intensity may be attributed to a lower excitation temper- 
ature (which is a reasonable assumption because CO has a much 
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Fig. 1. Left: Hubble Space Telescope (HST) image of the M 82 galaxy taken from http://hubblesite.org (Credit: NASA, ESA, and The 
Hubbl e Heritage Team (STScI/AURA)). The right panel shows the rotated HCO (2-1) interferometric map by Garc fa-Burillo et ail 
(2002) towards the central region of the galaxy. The dynamical centre is marked with a small black square. The north-east (NE) 
and south-west (SW) molecular lobes are shown with two dashed squares either side of the centre. The nominal position of our 
observations, at an offset (+13, +7.5), is marked with a cross. The beam sizes of the 2 mm survey for the lowest and highest 
frequencies are marked with two solid circles. The average beam size at 1.3 mm frequencies is marked with a dotted circle. 



lower dipole moment than most other molecules with detectable 
lines), and/or a lower beam filling factor (again something to be 
expected because the higher dipole species trace a denser gas 
component). 

Nevertheless, from these arguments it is plausible to assume 
that all the lines observed in our survey are optically thin or 
at most only moderately optically thick. The strongest lines in 
our survey are those of CS, C 2 H, and CH 3 CCH. The C 32 S/C 34 S 
intensity ratio of 13 (calculated using the (3 - 2) transitions 
detected in our survey) , is smaller than the isotopic ratio of 23 
derived for our Galaxy dWilson & Rood 1994), but there are rea- 
sons to believe that 34 S may have a higher abun dance in starburst 
galaxies than in our Milky Way (Ma rtin et al.ll2005l) . For C 2 H, 
a very str ingent lower limit t o its 13 C substitutions was deter- 
mined by Mar tin et~ai] ((2010), which suggests that the lines of 
CCH are also optically thin in M82. In addition, for the case of 
CH3CCH it is plausible that its lines are optically thin because 
this is a heavy molecule with a K-ladder structure where, hence, 
the population is spread over many energy levels. 



3.2. Source size 

On the b asis of M82 inte rferom etric observations of 12 CO 
and C ls O dWeifi et al.l20"oIah. HCO dGarcfa-Burillo etaD2002h . 
HOC + , and H 13 CO + dFuente et al.ll2008l) . we assumed an aver- 
age source size of 12" for all the molecular species detected in 
this survey. However, this is only a rough approximation because 
other species might exist in more clumpy or extended regions. 
To quantify the uncertainties in different source sizes for each 
molecule, we assumed a 9 S variation of 50%. We then recalcu- 
lated all the column densities and rotational temperatures for 9 S = 
6", 12", and 18". The emission size is very unlikely to be ei- 
ther smaller or larger than these values (although interferometric 
maps would be necessary to prove this). We found that the un- 
certainty in the source size has no strong influence on the derived 



N mo i and T mt parameters, limiting variations to maximal factors 
of 3.5 and 1.3, respectively. 



3.3. Different volumes 

Another uncertainty related to N mo \ and T Iot was introduced by 
our assumption that all the transitions of a certain molecule arise 
from the same gas volume. Low-/ transitions might arise from 
m ore extended region s than high-/ transitions. In a similar way 
to lBavet et aTld2009bl) . we selected CH 3 CCH (methyl acetylene) 
as a good example of a species for which we detected several 
transitions, because its lines are well spread over both the 2 mm 
and 1.3 mm bands. Thus, to quantify the contribution of this un- 
certainty, we plotted again the Boltzmann diagrams by succes- 
sively decreasing the source size by 10% for each higher tran- 
sition, i.e., CH3CCH (8 - 7) (which is the lowest transition we 
detected) had a S = 12", while CH3CCH (16 - 15) (the highest 
transition we have) had a 8 S — 7.2". When comparing the new 
Boltzmann diagram results with the ones for which we applied 
the same source size to all the lines, we found a difference of 
only ~ 14% in both iV mo i and T Iot . 



3.4. Rotational temperatures 

To accurately determine rotational temperatures from the 
Boltzmann plots, one needs at least two lines with sufficient 
energy differences between them. Otherwise, a large range of 
T Iot values satisfies the Boltzmann equation. This is indeed the 
case for some species for which we detected only one (/ - /') 
transition (e.g. C2H, H 13 CO + , or HCO). In these cases, we fixed 
r, ot = 20 + 10 K assuming that this is the average temperature 
for most of the observed molecules in the NE lobe of M 82 (see 
Table [TJ. This assumption might have an important impact on 
the derived column density for extreme temperatures well above 
or below this range (i.e. if T rot > 100 K or T rot < 5 K). 
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Fig. 2. M82 frequency survey coverage between 129.8 GHz and 175.0 GHz and between 241.0 GHz and 260.0 GHz. The spectra 
were smoothed to two channels, corresponding to a channel separation of ~28kms~ 1 in the 2 mm atmospheric window and to 
~18 kms -1 in the 1.3 mm atmospheric window. Species labeled with ;'- come from the image band. 
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Fig. 2. Continued 



4. Molecules in M82 

4. 1 . Detected and undetected molecules in the NE molecular 
lobe 

We detected a total of 69 molecular transitions and identified 1 8 
species, as well as 3 hydrogen recombination lines. While there 
are no newly detected extragalactic molecules, eight species 
were detected for the first time in M82. Furthermore, there 
are fourteen important molecular species which remained un- 
detected. Detected and undetected species are listed in Tables Q] 
and [2] respectively. The line ide ntifica tion was made usin g the 
CDMS (Miiller et afl l200lLl2005l) . JPL dPickett et al.ll 19981) . and 
Lovas (lLovasll 992) catalogues. To confirm whether a suspicious 
feature is either a weak line or a fake feature, we considered 
the following criteria: a) those features with line widths much 
narrower that the typical M82 line widths are unlikely to be 
weak lines in this galaxy; b) very complex molecules or very 
rare isotopologues (such as NH 2 CH 2 CH 2 OH, H 2 NCH 2 COOHII, 
35 C10N0 2 ... etc) are unlikely to be weak lines in M82; c) transi- 
tions with high energies of the lower state (usually > 100 cm"') 
are likely to be too weak to be detected in M82; d) lines with 
logio(integrated intensity) < -6 nm 2 MHz are likely to be too 
faint to be detected in M82. After ruling out cases of clear noise 
spikes, there were some remaining molecules with typical M82 



line widths, that were not too complex, of low E\ ow , and suf- 
ficiently strong. In each case, we examined all the transitions 
of the species that lie in the survey and checked whether other 
lines with similar spectroscopic parameters also seemed to be 
weak lines or not. This scrupulous procedure allowed us to dis- 
tinguish fake features and weak lines. We note that some of the 
spectra may contain ripples that were not completely removed 
because of the low degree of the baseline subtracted, which thus 
appear as weak lines. We also checked the image sideband for 
any possible lines. As a consequence of the high rejection, only 
H 2 CO (2 1 , 2 - 1 11) and HCO + (3-2) were detected from the image 
band at 133.2GHz and 259.4GHz, respectively. 



Gaussian profiles were fitted to all the detected lines. The 
parameters resulting from those fits are given in Table IA.1I In 
some cases, emission from the centre of the galaxy was detected 
by the telescope beam, displaying a bump at lower velocities 
of the lines. In these cases, we subtracted the bump before fit- 
ting the Gaussian profile. In other cases, two or more molecules 
were blended, thus we calculated the contribution of each line 
to the total profile and performed a synthetic fitting. The reduc- 
tion of the spectra and Gaussian profile fitting were performed 
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Table 1. Column densities, rotational temperatures, and fractional abundances of all the observed molecules in the NE lobe of M 82 



Molecule " 


N m0 ] (cm 


) 


1 rot 


(K) 


A'mol/A'c- 14 S 


C 2 H 


(1.3 + 


0.6) 


X 


10 15 


20.0 : 


t 10.0 


118.2 ± 60.7 


NO 


(1.2 ± 


0.1) 


X 


10 15 


8.8 : 


t4.4 


109.1 ± 10.7 


CH 3 CCH C 


(1.2± 


0.1) 


X 


10 15 


31.6 


±0.1 


109.1 ± 10.7 


CS(1) C 


(2.0 + 


0.1) 


X 


10 14 


5.9 : 


tO.l 


22.0 ± 1.7 


CS(2) C 


(4.2 ± 


0.1) 


X 


10 lj 


14.7 


+ 0.8 




CH 3 OH 


(1.2± 


1.1) 


X 


10 14 


26.2 : 


t 11.7 


10.9 + 11.7 


H 2 CO 


(1.1 ± 


0.2) 


X 


10 14 


25.5 


± 3.0 


10.0 ± 2.7 


H 2 S 


(6.1 ± 


4.5) 


X 


10 lj 


20.0: 


t 10.0 


5.5 ± 5.2 


C-C3H2 


(4.8 ± 


1.3) 


X 


10 lj 


10.1 


± 1.2 


4.4 ± 2.0 


HCO 


(3.9 ± 


2.8) 


X 


10 lj 


20.0 : 


1 10.0 


3.5 ± 3.5 


so f 


(3.7 ± 


0.5) 


X 


10 lj 


20.0 : 


1 10.0 


3.4 ± 1.2 


H 2 C» 


(3.1 ± 


2.5) 


X 


1 nl3 
10 


20.0 : 


t 10.0 


2.8 ± 3.2 


HC 3 N' 


(2.5 ± 


1.2) 


X 


10 


24.2 


± 3.9 


2.3 ± 1.9 


so 2 f 


(1.2 ± 


0.6) 


X 


10" 


20.0 : 


1 10.0 


1.1 ± 1.3 


NH 2 CN 


(1.2 ± 


1.5) 


X 


f All 

10" 


80.3 : 


t 52.9 


1.1 ± 2.2 


c 34 s 


(1.1 ± 


0.7) 


X 


10 13 


20.0: 


t 10.0 


1.0 ± 1.0 


CH3CN 


(6.4 ± 


1.2) 


X 


10 12 


32.9 


± 1.6 


0.6 ± 0.8 




(1.7 ± 


1.6) 


X 


10 12 


20.0: 


t 10.0 


n _t_ n o 
u.z + u.y 


H 13 CO + 


(1.3 ± 


0.8) 


X 


10 12 


20.0: 


t 10.0 


0.1+0.8 


c i* d 


(3.5 ± 


1.7) 


X 


10 16 


20.0: 


t 10.0 


3181.8 + 1700.7 


CN" 


(2.0 ± 


0.5) 


X 


10 14 






18.2 + 5.7 


HCN" 


(4.0 ± 


0.5) 


X 


10 13 






3.6+1.2 


HNC / 


2.2 x 1 


o 13 


10.0 


2.0 


CO +s 


(6.6 ± 


3.4) 


X 


10 12 


20.0: 


t 10.0 


0.6+1.0 


HOC + h 


(3.9 ± 


1.9) 


X 


10 12 


20.0: 


t 10.0 


0.4 + 0.9 



Molecules in boldface are new detections in M 82. 

Some molecules do not have a wide enough dynamic range in energy to calculate the rotational temperature with precision. For them, we 
assumed T Iol = 20 ± 10 K (see Sect. l3.4l an d AppendixlAl for details). 
Some extra detections observed by I Aladro et alj d201 lah were used. 

CS shows two components with different column densities and rotational values corresponding to the two fits that result from its Boltzmann 
diagram. The fractional abundance value is in this case [A^^s^j + ^ cs(2)]/7Vr 34 s- 

The mo lecules shown at the botto m of the table were observed by: d Martin et al. (2009a); e Fuente et al. (2005) (averaged value for the galaxy 
disk); ^ iHuttemeister et all dl995h (averaged value for the n uclear zone); s Fuente et alJd2006t) : * Aladro et al. in prep. 



Only upper limits were reported until now: SO, and S0 2 by Petuchowski & Bennettl d!992l) . and H I3 CN by I Wild et al] d!990h . 



using the CLASS and MASSA0 software packages. Each case 
is described in Appendix lAl 

From the LTE analysis, we found that C2H is the most 
abundant molecule in our M 82 survey, followed by NO, which 
is detected for the first time in this galaxy. The 13 C isotopo- 
logues H 13 CO + and H 13 CN have the lowest column densities 
in the M82 NE molecular lobe (see Table[TJi. We augmented our 
CH^CCH, CS and H C^ N data with other det ections observed by 
lBavetetal.l(l2009bh and lAladro et al.l(l201 lal) in frequencies that 
lie outside our frequency coverage. For the whole list of detected 
molecules, CS is the only one that has two clear different rota- 
tional temperatures (of ~6 K and ~15 K). In addition to CS, NO 
also has a very low rotational temperature, being <10K. The 
molecule NH2CN could have a high rotational temperature (of 
~80 K), th at is similar to what is observed in the starburst galaxy 
NGC 253 dMartfn et al.ll2006bh . 

4.2. Fractional abundances 

We calculated the fractional abundances of each species with re- 
spect to C 34 S (Table [TJ, and the upper limits to the abundances 
for the undetected species (Table [2]). We used the C 34 S (3 - 2) 



2 CLASS http://www.iram.fr/IRAMFR/GILDAS 

3 MASSA http : 1 1 damir.iem.csic.es I mediawiki — 
1.12.0/ index, ph p /MAS SAM ser' s Manual 



line for this purpose because it is supposed to be significantly 
more optically thin (r << 1) than the main isotopologue, which 
in the worst case could be m oderately saturated, with r < 2 
dMauersberger & Henkellll989h . Nevertheless, another abundant 
optically thin gas tracer, such as C I8 0, could also have been 
used. We note, however, that C 18 arises from more extended 
regions, with critical densities of at least one order of magnitude 
lower than those traced by C 34 S. Nonetheless, C ls O is some- 
times considered the most effective way to derive H2 column 
densities. The problem here is that the 16 0/ 18 ratio has not yet 
been tightly con strained, having on ly a lower limit of 16 0/ 18 
> 350 for M 82 dMartfn et al.ll2010h . As a consequence, the H 2 
column density cannot be accurately determined. The same ar- 
gument can be applied to C 34 S, because the s ulphur ratio 32 S/ 34 S 
is not well-known either dMartfn etalJ l2010). Nevertheless, we 
prefer to refer our fractional abundances to the column densities 
of a dense gas tracer such as C 34 S, because CS and its 34 S iso- 
topologue have similar column densities in M82 and NGC 253 
(see Sect.|5]for more details). 

4.3. A picture of the nuclear region of M82 

To have a more complete picture of the chemistry in the cen- 
tral region of M82, we compiled from the literature the results 
of other molecular detections observed towards the central hot- 
spot and the surrounding molecular ring. As mentioned in Sect|2] 
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Table 2. 3cr upper limits to the column densities and fractional 
abundances of some important undetected molecules. 



Molecule 


N (cirr 2 ) 


N/N c m s 


ocs 


< 7.4 x 10 13 


< 67.0 x 10-' 


HNCO 


< 2.5 x 10 13 


< 23.0 x 10-' 


C 2 S 


< 1.0 x 10 13 


< 9.0 x 10-' 


c-C 3 H 


< 9.7 x 10 12 


< 8.8 x 10-' 


C 2 D 


< 9.2 x 10 12 


< 8.4 x 10-' 


CH 2 NH 


< 8.4 x 10 


< 7.6 x 10~' 


13 CS 


< 4.8 x 10 12 


< 4.4 x 10-' 


HOCO+ 


< 4.2 x 10 12 


< 3.8 x 10-' 


NS 


< 2.3 x 10 12 


< 2.1 x 10-' 


HN 13 C 


< 2.3 x 10 12 


< 2.1 x 10" 1 


SiO 


< 1.6 x 10 12 


< 1.4 x 10-' 


DCN 


< 1.2 x 10 12 


< 1.1 x 10-' 


DNC 


<7.0x 10" 


< 3.4 x 10~ 2 


N 2 D + 


< 5.3 x 10" 


< 4.8 x 10" 2 



To calculate the upper limits to the column densities, we assumed 
r rot = 20 K and Av = 100 km s _I for all the molecules. 
Fractional abundances are compared to those derived for NGC 253 
in Fig. [5] 



our observations were carried out toward s the NE molecular 
compl ex, because the HCO observations of iGarcfa-Burillo et al.l 
( 2002) show that the photodissociating radiation is stronger there 
than in other parts of the nucleus. The emissi on in the NE molec- 
ular lo be is also stronger when observing CS dBavet et al.l20 08b, 
2009b), which is claimed to be a good tracer of dense gas linked 
to the presence of UV radiation. However, other species such 
as CO dMao et al.ll2000l) H 7 CO (iHiittemeister et al.ll 19971) . and 



HNC (Hiittemeister et alJ 19951) s eem to be more abundant i n the 
SW lobe, while N 2 H + and C 3 H 2 dMauersberger et al.ll99lh may 
peak at the very centre. In particular, it was found that a high 
abundance of CO in combination with a high gas density may 
lead to a rapid destruction of N 2 H + in the molecular ring via the 
reaction N 2 H + + CO -> N 2 + HCO + dMauersberger et allll 99 U 
J0rgensenetal . 200|). 

Each molecule may arise from gas that can be at different 
temperatures. There is a notable difference between the tem- 
perature derived from two of the most important thermome- 
ters, H 2 CO and NH 3 . iMiihle et all d2007l) carried out a multi- 
transition analysis of formaldehyde in the NE and SW molecu- 
lar lobes. They derived densities of 7 x 10 3 cirT 3 and high ki- 
netic temperatures around 200 K, and found that H 2 CO seems to 
be co-spatial with CO. The relatively low densities obtained for 
both species are qu ite similar. However, the study performed by 
IWeifi etalJ d2001b) towards the SW molecular lobe using ammo- 
nia indicates that the NH3 emitting gas is rather cold, with kinetic 
temperatures of ~ 60 K. The fractional abundance of NH3 is ex- 
tremely low in M82, because this molecule is easily destroyed by 
the UV radiation. If we were to assume that the ammonia prop- 
erties in the NE lobe are similar to those in the SW, it could be 
said that most of the observed molecular gas can be more accu- 
rately represented by H 2 CO rather by NH3 . On the other hand, 
the low density determined from the H 2 CO in M82 could lead 
to highly subthermal excitation. Furthermore, we expect that at 
least part of the detected molecular emission arises from well 
quite-shielded molecular regions, which also produce NH3 emis- 
sion, and should be relatively cool. Thus, the kinetic temperature 
should be considered in each single case to see whether the prop- 



erties and spatial distribution of a given molecule are more sim- 
ilar to those of formaldehyde or ammonia. 

5. Comparison of the M82 and NGC 253 chemistries 

The galaxies M82 and NGC 253 share several characteristics 
that have ensured that they are among the most studied of ex- 
tragalactic objects. Apart from their luminosities (Lf, GC :x ' - 



2.1 x 10 10 LM.Ideyaucouleurset al.lll99"Tt Z™ 82 = 3 x 10 10 L Q , 

lie at 



TelescoJ 



iRice et alJ [1988T 
tance (D NG c253 = 3.9Mpc, D M 82 
1994), so their morphology, physical processes 



they lie at a simi lar dis- 
; 3.6 Mpc, iFreedman et al] 
and chem- 



istry can be studied in detail with the same spatial resolution. 
Furthermore, both nuclear regions host very active star-forming 
regions (SFIW953 ~3.6M yr"\ SFR M8 2 ~ 9M yr-', 
Strickland et al. 2004). 



Martin et al. (2006b) carried out a frequency survey towards 
the central region of NGC 253 covering the same frequency 
range in the 2 mm atmospheric window as the M82 survey 
presented here. Up to 25 species were identified in NGC 253. 
The combination of both surveys allows us, for the first time, 
to perform a detailed and complete comparison of the chemi- 
cal composition in these prototypical starbursts. Figure [3] shows 
the comparison of the 2 mm scans towards both galaxies. In 
Fig.|4] we compare the fractional abundances of the 18 species 
detected in both line surveys. We also included other impor- 
tant m olecules f rom the literature, such as C 18 Q (Martinet all 



2006), 



and 



l2009allbl). CO + dMartfn et al] l2009bl; iFuente et al 
HOC + ((Martin et alj| |2009bl Aladro et al. in prep.) since they do 
not have transitions in the covered frequency range. Moreover, 
we also show in Fig.[5]a comparison of the fractional abundances 
of 12 other important molecules that were not detected in our 
M 82 survey, but towards NGC 253. Therefore, we plot 3<x upper 
limits to their abundances. In addition, Fig.|6]compares the rota- 
tional temperatures of several molecules in M82 and NGC 253. 

Next, we discuss in detail all the species, emphasizing those 
that allow us more clearly differentiate between the M 82 and 
NGC 253 chemistries. 

- C 2 H. The ethynyl radical may be strongly related to the 
Ci abundance and its fo rmation can be exp lained in terms 
of gas-phase chemistry dTurner et al.l lT999). The main for- 
mation path of C 2 H is via electron rec ombination re- 
actions with the C ?H:f and C?Uj ions dWatsonl Il974t 
IWoottenetal.lll980t iMul&McGowanl Il980t) . It has been 
claimed to be a good photon-dominat e d region tracer 



(Huggins et al 



19841: iHenkel et al] 119881: iMeier & Turned 
120051; IFuente et al.ll2005t iBavet et al]|2009bl) . The molecule 
C 2 H is a factor ~ 3 more abundant in M 82 than in NGC 253. 
Furthermore, it is the brightest spectral feature in both galaxy 
surveys. 

Although this molecule was detected in a large number of 
Galactic sources, it went almost unnoticed in external galax- 
ies, except for the M82 and NGC 253 detections (see also 
Martin et al. 2010), up to now it has only been observed to- 
wards one additional galaxy, namely IC 342. There it was 
found in th e central 50- lOOpc re gion, where radiation fields 
are strong dMeier & Turneril2005h . 

CS. Carbon monosu lfide is one of the most well-studied 
dense gas tracers (Ba vet et ail l2008al) in extragalactic 



sourc e s, and in particular starburst gala x ies ( Henkel & Ballv 
1985; Mauersberger & Henkel 1989; Mauersberger et al. 
ri989t IBavet et al.ll2008bl l2009bt FAladro et al.l l2011al) . The 
main chemical reactions that lead to the CS formation are 
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Fig. 3. Comparison of the M82 and NGC253 surveys at 2 mm wavelengths. Both spectra have been smoothed to a velocity res- 
olution of ~28kms~ 1 . At first glance, NGC253 shows a greater chemical complexity than M82, with more lines and molecular 
species. 



detailed in Pineau des Forets et al. I (119861) and DPrdla et al. 
(fl989h . A high CS abundance is predicted for slightly 
shielded regions dSternberg & Dalgarnolll995l) . The CS col- 
umn densities and fractional abundances with respect to 
C 34 S are almost the same in both M 82 and NGC 253, which 
is unsurprising if CS saturation is insignificant (or compara- 
ble), and the 32 S/ 34 S isotope ratios (both sulphur nuclei are 
for med by oxygen burning in m assi ve stars) are similar (see 
alsolReauena-T orres et alJl2006l and lMartfn et alJ l2009a). 
It can be seen in Fig. [6] that CS has up to three r iot compo- 
nents for NGC 253, while M 82 only shows two components 
(Fig. E2 Table [1). The highest r rot value of NGC 253 is 
well above the highest one in M 82. This could be explained 
if the gas in NGC 253 had a higher kinetic temperature. It 
seems that in these two galaxies, the spatial distribution of 
CS is similar to that of a mmonia, or at least, i s clumpier than 
those of H 2 CO and CO d Aladro et al.1 1201 lab . If so, the ki- 
netic temperature in M 82 would be around 60 K, while that 
of NGC 253 would be between 100 and 120 K. In that case, 
NGC 253 does not necessarily need to have higher H2 densi- 
ties than M 82, as seen by the LVG modellin g of the CS emis- 
sion towards these two galaxies bv lAladro et al.l d201 lab . 
As mentioned in Sect. 14.21 C 34 S is strong enough to be 
clearly detected in both galaxies and is observed to be sim- 
ilarly abundant. On the other hand, although 13 CS was de- 
tected in NGC 253, it is fairly faint in M 82, where we only 
have an upper limit. Finally, the double isotopologue 13 C 34 S 
is too faint to be detected in both galaxies. A comparison 
of the CS (3 - 2) and C 34 S (3 - 2) integrated intensities indi- 
cates a sulphur isotopic ratio of 32 S / 34 S ~ 12. This value may 
be hi gher, since CS could be moder ately optically thick in 
M 82 (Mauersb erger & Henkellll989h . Using the upper limit 



to 13 CS (see Table [2}, we also calculated a lower limit to the 
12 C/ 13 C carbon isotopic ratio. We obtained 12 C / 13 C > 50. 
This h mit is lower than the one obtained by Marti n et~a"D 
(2010) using the molecule CCH and its carbon isotopologues 
( 12 C/ 13 C > 138). 

C-C3H2. As for CS, cyclopropenylidene is also formed by 
means of gas phase ion-molecule reactions, from the dis- 
sociative recombination of the stable cyclopropenyl cation, 



C3HI and its cyclic and linear isomers, C-C3H 3 and I-C3H3 



(Thad deus et alJI 19851 lA"dams & Smithj|1987l) . Thus, the de- 
struction of some species by UV fields and the subsequent 
creation of ions would favour an enhancement of the C-C3H2 
abundance in M82. Our M82 data indeed indicate a frac- 
tional abundance four times higher than in NGC 253. The 
rotational temperature derived from the Boltzmann diagram 
of M 82 is 10.1 + 1.2 K (similar to the temperatures obtained 
in NGC 253). This seems to indicate that cyclopropenylidene 
arises from the envelopes of the molecular clouds. 
H 13 CO + , H I3 CN. These two species are the only I3 C iso- 
topologues detected in our survey, and have the lowest abun- 
dances among all the molecules detected in this study. The 
ratio of their main carbon isotopologues, HCO + /HCN, has 
been studied to differentiate between the physical processes 
occurring in starburst galaxies and AGN envir onments, and 
has been found to be lower in the latter dKohno et al.l 
|200H llmanishietall 12001 120091 iKrips et alJ l2008h . The 
HCO + /HCN ratio also reflects how evolved a starburst re- 
gions is, being higher in evolved starbursts (where there 
are strong ionization effects) than in young starbursts (that 
are more dominated by shocks). This ratio can also be 
calculated using the 13 C isotopologues, with the advan- 
tage that the 13 C lines are likely to be optically thinner, 
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Fig. 4. Comparison of the M82 and NGC253 fractional abundances of different molecules with r espect to the C 34 S (3- 2) abun- 
dances. All the fractional ab undances have been calculated from the millimeter freque ncy surveys oflMartfn et al.l d2006bl) and this 
paper e xcept for: (a) C 18 Q: | Martm et alld2009al): (b) C 18 Q HCO, HOC + , and H 13 CO + : lMartm et al.l(l2009bl) : (c) CTUCCH. CS and 
HC.N: lAladro et all (1201 lah : (d) CO + : lFuente et al.l (120061: (e) HOC + : Aladro et al. in prep. 



thus avoiding line opacity effects. From our data, we ob- 
tain H 13 CO + (2 - 1)/H 13 CN(2 - 1) = 1.9. This is signif- 
icantly higher than HCO + (l - 0)/HCN(l - 0) = 1.4 and 
HCO+(3 - 2)/HCN(3 - 2) = 1.2, derived by iKrips et all 
(2008) using the main isotopologues for the same region 
in M 82. However, the difference in the ratios HCO + /HCN 
and H 13 CO + /H 13 CN may also be due to opacity effects since 
the HCN line could be less saturated (because of the hyper- 
fine structure) than the HCO + lines. On the other hand, our 
H 13 CO + (2 - 1)/H 13 CN(2 - 1) ratio for M 82 is a factor of 
~ 5 higher tha n the one found fo r the same transitions in 
NGC 253 (OA lMartfn et al.ll2006bh . This difference between 
M 82 and NGC 253 is consistent with the idea that evolved 
starbursts have higher HCO + / HCN ratios. 
HCO, HOC + , CO + . These molecules are claimed to 
be good PDR tracers, with their abundances being en- 
hanced in zones where the ch emistry is strongly influ- 
enced by UV and X-ray radiation dG arcia-Bu rillo et al.l2002b 



Savage & Ziurvs 2004; Fuent e et al.ll200 5. 2008; Gerin et alj 
2009; Mar tin et al.ll2009bl) . Thus, these species are formed 
through ion-molecule reactions in gas phase processes, al- 
though the chemical path of HCO is l ess clear and could 
also involve neutral-neutral reactions dWatsonetal.il 1975b 
iHollis & Churchwelllll983HSpaans & Meiierinkl 120071) . The 
HCO-to-H 13 CO + ratio has been found to be a good indica- 
tor of UV fields in the before mentioned zones. From our 
data, we find that A^hcop-i) / N-Hnr. n+o.-n _ 30. This va lue is 
similar to those found in N GC 253 (Ma rtin et al.l2009bl) . and 
Orion (ISchilke et al.ll200lT) . Ho wever, this ratio is a factor of 
ten higher than the one found bv lGarcfa-Burillo et al.l (2002) 
for the whole disk of M82 (A^HCOd-O) /-Wh°co+(i-0) ~ 3.6). 
On the other hand, the ratios of CO + and HOC + with re- 
spect to HCO + calcu lated towards both galaxies are similar 
dMartfn et al.ll2009bl) . 

To complete our sample of molecules, we included the 
HOC + (l - 0) and CO + (5/2 2 - 3/2j) lines from lMartfn eTaTI 



d2009bl) . iFuente et all (120061) . and Aladro et al. in prepara- 
tion (see Table [TJ and Fig. |4]for details). As expected, both 
molecules are more abundant in M 82 than in NGC 253. This 
is especially true for CO + , whose fractional abundance is al- 
most four times higher. The large amount of this short-lived 
reactive ion is clearly indicative of the strong PD Rs towards 
the NE molecular lobe of M 82, as shown by IFuente et al.l 
(2006). 

- CH3OH and HC3N. Methanol and cyanoacetylene can 
be efficiently formed on grain mantles rather than 
by means of gas -phase reactions (Harta uist et al.l 11995b 
Hida kaeTall 12 004). and ar e claim ed to be good tracers of 
dense gas ( Rodrigue z-Franco et all 1 998b iBavet et alj2 008a; 
lAladro et al.ll2.01 lah . Methanol can also trace shocks when it 
is released from the grain mantles to the gas phase through 
grain sputtering in shocks. Both species are supposed to be 
easily dissociated in the external and less dense layers of the 
molecular clouds. At the same time, their abundances in- 
crease in the molecular cloud cores, where the gas is well- 
protected from the UV radiation. Both CH 3 OH and HC 3 N 
have fractional abundances 2-3 times higher in NGC 253 
than in M 82. This is consistent with the results of lSnell et all 
(120111) . 

The temperature and colunm density that we obtained 
for methano l in M 82 are similar to the ones found by 
Martin et al. (2006a) for the densest methanol component 
detected in the NE molecular lobe of M 82. However, we 
do not find evidenc e of the colde r and le ss dense gas com- 
ponent reported by Ma rtin et al.1 (l2006al) (with Nqh^oh — 
2.7 x 10 13 crrT 2 and T rot = 5.0 + 0.1 K), since this compo- 
nent is traced by a lower excitation transition at 3 mm wave- 
lengths 

Methanol generally arises from cold gas in galaxies 
dMartfn et al.ll2l)06allbT) . which is in good agreement with the 
NGC 253 rotational temperature seen in Fig. [6] However, the 
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Fig. 5. Fractional abundance limits of twelve undetected 
molecules in our M 82 survey compared to the N GC 253 abun- 
dances taken from the frequency survey done by Martin et al. 
(2006b). HNCO abundances were taken from iMartfn et al.l 
h009al) . 



Trot in M 82 is warmer (~ 26 K), leading to a difference in 
rotational temperatures between both galaxies of ~ 15 K. 
On the other hand, cyanoacetylene u sually traces very ex- 
cited and warm gas (Aal to et aljlfool . While this molecule 
identifies a constant temperature distribution in M82 (T TOt ~ 
24 K), it shows up to three density components with different 
rotational temperatures in NGC 253 (the component belong- 
ing to the higher excitation transitions has a rotational tem- 
perature of ~ 73 K). This is due to the different density struc- 
ture of the molecular clouds in both galaxies (Aladr o et al.l 
1201 lah . The molecule HC3N indeed displays the largest dif- 
ference in T rot between M 82 and NGC 253 (AT mt ~ 50 K, 
for the warmest NGC 253 component, Fig. g). That HC 3 N 
has a very warm T rot in NGC 253 (73 K) seems to be related 
to the higher kinetic temperature of this galaxy, since the 
averaged nH, densiti es are lower in NGC 253 than in M82 
dAladroetal.ll2011al) . 

The [C-C3H2/HC3N] ratio appears to in dicate the evolution - 
ary stage of the bursts of star formation dFuente et al.l F2005). 
Since C-C3H2 is enhanced in PDRs, and HC3N is instead 
destroyed in these regions, this ratio should be higher in 
M 82 than in NGC 2 53. Comparing our results with those of 
Marti neTaT] d2006bl) . we find a difference amounting to one 
order of magnitude with N c -c 3 n 2 1 Nhc 3 H - 2.0 and 0.2 for 
M 82 and NGC 253, respectively. 
- CH3CN. Methyl cyanide can be for med via gas-phase, 
where CHt and HCN a re precursors dMauersberger et all 
fl99lt iMiliar et al.l \l99% . These molecules CH 3 CN and 
CH3CCH are goo d tracers of dense gas that h a ve been used 
as the rmometers (IChurchwell & Hollisl 119831: iGusten etall 
fl98l . The difference in their dipole moments (0.75 Debye 
for CH3CCH and 3.9 Debye for CH 3 CN) can lead to dif- 
ferences in the r iot , because molecules with low dipole mo- 
ment are expected to trace temperatures closer to the kinetic 
temperature. However, we find similar rotational tempera- 



tures for both molecules in M82 (T rot = 32.9 + 1.6K for 
CH 3 CN and r rot = 3 1 .6 + 0. 1 K for CH3CCH). The expected 
behaviour is observed in NGC 253, where CH3CN shows 
a lower T m than CH3CCH (9.6 ± 0.7 K and 63.0 + 20.0 K 
respectively. IMartfn et al.ll2"006bl) . Thus, the relationship be- 
tween the dipole moments and r rot is unclear. In any case, 
CH3CN arises from a more clumpy gas, possibly from the 
cores of the molecular clouds, while CH3CCH, more abun- 
dant and extended, a rises from intermediate density regions 
( Alad roet al.l 1201 Tab . Furthermore, CH3CN seems to be 
tracing gas with different properties in both galaxies: the 
molecular cloud cores of M 82 traced by this species seem 
to be warmer than those of NGC 253 (Fig.|6]l. 
Sulphur-bearing species. In general, sulphur is believed to 
freeze onto grain mantles, and experience a later evapora- 
tion by shocks that leads to the subsequent formation of sev- 
eral sulphur species, such as PCS, SO, SO2, Hg S, or H2CS 
dScalo & Slavskvl Il980t IMiliar & Herbs! fl990h . However, 
sulphur might frequently alternate between a grain-phase 
and a gas-phase because of the volatile and fr agile nature 
of some of the molecula r compounds it forms (Minh et al. 
1199 It lOmont et al.lll993l) . The aforementioned species are 
not only shock tracers but also tracers of high-mass star for- 
mation in starburst ga l axies s uch as M 82, as proved by the 
models of iBavet et all (l2008al) . 

We have found that SO, SO2, and H2CS have similar abun- 
dances in M82 and NGC 253. This indicates that exist 
shocks in the NE molecular lobe of the galaxy. We know 
that M 82 is not a pure-PDR but it also contains signific ant 
amounts of dense UV-protected gas (Martin et al. 2006a) af- 
fected by shocks, as reported by Garcfa-Burillo et al. (2001) 
using SiO observations. The low SiO abundance in M82 
(in comparison with that of NGC 2 5 3) was firstly measure d 
by iMauersberg er & Henkell (1 19911) . ISage & Ziurvsl dl995h . 
and iGarcia-Burillo et al.l (120001) Similarly, a lower PCS 
abundance was reported in M82 ((Mauersberger et al.lll995h . 
which we confirm here by its non-detection in our line sur- 
vey. 

H2CO. Formaldehyde can be efficiently formed on the 
surface of icy grain s via the hydrogen afion of CP 
dCharnlev et al.1 119971: [Watanabe & Kouchil |2002|). It is 



2008a). 



claimed to be a tracer of dense gas (Bav et et al 
In those occasions when it is released into the gas- 
phase by shocks (thus also behaving as a shock tracer), 
it has abundances simi lar to those of other molecules 
such as methanol or SiP (Bachille r & Perez Gutierrezll 19971: 



IRodriguez-Fernandez et aL 20101) . The H?CG fractional 
abundances (relative to C 34 S) of M 82 and NGC 253 are 
similar, although that of M 82 might be the consequence of 
photo-evaporation rather than shocks. 
The rotational temperatures of H2CO are similar in NGC 253 
and M 82 (~ 26 and 31 K respectively, see Fig.[6]l. 
HNCO. Chemical models include dust-grain surface re- 
actions and gas-phase reactions for explaining the forma- 
tion of isocyanic acid (HN CO), because gas-phase reactions 
alone seem to be inefficient (ICaselli et al.l 19931: iGarrod etall 
120081: IMartfn et al.l l2009bl) . The molecule HNCO is sup- 
posed to be released through grain mantle evaporation by 
shock waves and/or by neutral- neutral reactions in the gas 
phase involving CN and O2 (Rodri guez-Fernandez et al.l 
2010) .This molecule has emerged as one of the most likely 
promoting star formation in the mole cular clouds of both 
the Milky Way and external gala xies dMartfn et al.l l2009a: 
IRodriguez-Fern andez et al. 2010). Its abundance varies by 
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Fig. 6. Comparison of the M 82 and NGC 253 rotational temper- 
atures for several molecules. 



nearly two orders of magnitude a mong galaxies dom inated 
by different physical processes dMartfn et alj 12009a) . Our 
data support this idea, because this molecule is not de- 
tected in our survey, while there is a clear detection in 
NGC 253. The molecule HN CO has never been d etected to- 
wards M 82 (see also iNguven-O-Rieu et alJI 199 lit , presum- 
ably b ecause this molecul e is very sensitive to photodissoci- 
ation (iMartfn et al.ll2009bl) . 

- C-C3H. Although the production mechanism of cyclopropy- 
nylidyne is not well-known, it has been proposed that it 
has the same ion precursor as C3H2, namely C3H3 , in dar k 
clouds dAdams & Smithlll987t iMangum & Woottenl fl990). 
In PDRs, the C3H2 abundance is seen to be higher than that 
of C-C3H. As an example, in the Horsehead nebula the C-C3H 
fractional abundance is one ord er of magnitude lower than 
the former dTevssier et al.ll2004l) . Though C-C3H could be re- 
lated to local PDRs, in M 82 it is too faint to be detected. We 
find the ratios C 2 H/c-C 3 H > 100 and c-C 3 H 2 /c-C 3 H > 5. The 
upper limit to the abundance derived from our non-detection 
is still similar to the value derived from NGC 253. 

- NS. We describe this sulphur-bearing molecule in more de- 
tail because it shows the largest fractional abundance con- 
trast between M 82 and NGC 253 after HNCO and CH3CCH. 
This mole cule was det e cted f or the first time outside our 
Galaxy by [Martin eTail d2005l) in NGC 253, and since then 
no other extragalactic detection has been reported. It was 
claimed that the NS chemistry is driven by the large-scale 
shocks that dominate the NGC 253 galaxy nucleus. On the 
other hand, the NS / CS ratio was pr oposed as an i ndicator of 
the presence of shocks in hot cores dViti et al]l200ll) . We ob- 
tain NS/CS < 0.01 and NS / CS = 0. 1 for M 82 and NG C 253 
respectively, using the data from Marti n et al ] d2006b). This 
difference of at least one order of magnitude indicates that 
UV fields affect the chemical abundance of NS, and indicates 
that this ratio is also a good indicator of shocks in galaxies. 



HOCO + . Protonated carbon dioxide (HOCO + ) is thought to 
be for med in the gas phase from desorbed CO2 dTurner et alJ 
[1991 . In our Galaxy, it was observed towards translucent 
clouds and the Gala ctic centre, being especially abundant i n 
SgrA and SgrB2 dMinh et all [19881: iDeguchi et1il] [2006). 
The enhancement of this molecule in these regions was 
linked to shock activities that release the CO2 frozen onto 
grains into the gas-phase. Outside the Milky Way, NGC 253 
is the only galaxy towards which this molecule has so far 
been detected. Our upper limit in M 82 shows that it is at 
least 2.5 times less abundant than in NGC 253. This under- 
abundance might be due to a lower influence of shocks in 
M82. 

CH3CCH. After HNCO, propyne is the molecule that has 
the largest difference between the abundances of NGC 253 
and M 82, which are eight times higher in the latter galaxy. 
This difference in the abundances was also found by 
ISnell et all d2011l) , who detected CH3CCH in M 82 but not 
NGC 253. The column density that we derive in M 82 (1 .2 x 
10 15 crrT 2 ) is the highest ever measured for an extragalactic 
source. This results in a fractional abundance CH3CCH/H2 = 
1.5 x 10~ 8 , that is similar to the abun dance found towards 
Orion (between ~ IP" 8 and ~ 10~ 9 . IChurchwell & Hollisl 
119831: ISutton et al] 119851) . The most plausible explanation 
of the high abundance of this molecule in M82 is that it 
was formed in the gas phase, in which the destruction of 
propyne by UV fields could be balanced by efficient for- 
matio n by means of ion-molecule or neutral-neutr al reac- 
tions dMillar & Freeman|[l984t iBisschop et al.ll2007l) . In ad- 
dition, propyne is observed in dense regions that are well- 
shielded from the dissociating radiation. Fig.[6]shows that the 
CH3CCH rotational temperature is more than 10 K higher 
in NGC 253 than in M82. This is due to a possible highe r 
kinetic temperature in NGC 253 (Maue rsberger et a"lH 2003). 
as well as a higher H2 densit y in the molecular clouds of 
NGC 253 traced by CH3CCH jAladro et alj201 lal) . 
NO. Nitric oxide is detected for the first time in M82. It 
is the main precursor of the N/O chemical network, and 
its abundance is directly related to the nitrogen and oxy- 
gen abundances by the reactions NH + O — > NO + H 
dHalfen et al.|[200Th and N + OH -> NO + H dLique et alj 
120091: iBergeatet aUl2011h . 

We note that the nitric oxide abundances are almost the same 
in M82 and NGC 253. These values indicate that it is not 
strongly affected by UV fields or shocks, as in the case of 
CS. On the other hand, NO is likely to be widespread across 
the nuclear regions. The NO dipole moment is quite similar 
to that of CO and it is therefore likely to trace the overall 
gas. The low rotational temperatures of NO in both galaxies 
(T mt = 8.8 ± 4.4 K for M 82 and 6 + 2 K for NGC 253) also 
indicate that it possibly originates in the external and colder 
layers of the molecular clouds. 

NH2CN. The formation process of cyanamide is not well- 
known, since this molecule was observed only towards a few 
sources. Its column density in M82 is con sistent with that 
found in SgrB 2 bvlCummins et alj d!986l) (A^nh,cn = 2 x 



10 13 cnT 2 ) and lTurneret all d 19751) (N NHiC n ~ 10 14 cnr 2 ). 
However, cyanamide is so faint in M 82 that we only have 
tentative detections. As a consequence, its rotational temper- 
ature and column density have very large errors, so these re- 
sults should be taken with caution. The molecule NH2CN 
was not detecte d in other local massive star-forming re gions 
such as Orion dSutton et alJll985h . iTurner et all dl975h con- 
cluded that this molecule probably arises from high density 
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regions with n> 10 5 cm~ 3 . In exter nal galaxies, cyan amide 
was also detected towards NGC 253 (IMartfn et al.l2006bl) . In 
contrast to the SgrB2 molecular cloud (OH position), where 
the excitation temperature of this molecule was found t o 
be relatively low (T rot = 15 + 3 K. ICummins et alj Fl986). 
cyanamide shows a very high rotational temperature in 
NGC 253 (T rot = 67 + 13 K). The latter v alue is closer to 
that f ound towards the SgrB2 (M) hot core (Nummelin et al. 
2000). This suggests that this molecule might arise from the 
very dense hot cores in starburst regions. However, our data 
can neither support, nor refute this possibility. 
CH 2 NH. Methyleneimine (or CH 2 NH) is possibly formed 
by neutral-neutral gas-phase reactions, although grain- 
chemistry ca nnot be ruled out. This species is indeed poorly 
understood dTurner et al.l 1 1999b . CH 2 NH was observed to 
be more abunda nt in Galactic star- forming regions than 
in dark clouds (iDickens et al.l 1 19971) . This indicates that 
this molecule is enhanced in warm environments, so it 
might be considered a good tracer of high temperatures. 
Outside the Mil ky Way, methylen eimine was detected to- 
wards N GC 253 dMartfn et al.l2006bl) and the ULIRG galaxy 
Arp220 dSalter et al.ll2008l) . Our non-detection in M 82 indi- 
cates that its abundance is at least a factor of 2.5 lower than 
in NGC 253. 

Deuterated molecules. No deuterated molecules have been 
detected in M 82. In accordance with the lBavet et al.l d2010l) 
models, at least D 2 CO, DCN, C 2 D, and HDCS should be de- 
tectable in both M82 and NGC 253. However only PNC, 
and N 2 D + are tentatively detected in NGC 253 (Marti n et al 
2006b). Using the HNC data obtained by Hiittemeist er et al 
(1995) for M82 and our upper limits for DNC, we obtain 
D/ H < 5 x 10~ 4 Simila rly, using the N 2 H + data obtained 
bv ISage & Ziurvsl d 1995b we obtain D/H < 5 x 10~ 2 . The 
limit derived with DNC is con sistent with the interstellar 
medium value of (0.5-2)x 10~ 5 ([Turner & Zuckermanl 19781; 
Walmslev et al.lll987t. iMauersberg er et alJll988HParise et al.1 
20091) . 



- Radio recombination lines indicate the presence of ionized 
gas and can be used as an indicator of the star formation 
activity. There is a direct relationship between the integrated 
emission of the H II regions and the O B star formation rate 
(ICohenlll976llKennicutt & KentHl983l) . Comparing our data 
and that of IMartfn et a0 (2006b), we find that the integrated 
emission of the Ha lines is more than 1.5 times higher in 
NGC 253 than in M 82. This indicates that there are more hot 
stars in the nucleus of NGC 253. It could, however, also be 
related to the evolutionary stage of the star forming regions 
at the centres of both galaxies, where young stars emitting 
radiation are more common in a early-to-intermediate stage 
of bursts of star formation than in old starbursts. 

6. Conclusions: A chemical scenario of the M82 
starburst 

At first glance, Figs. |4] and [5] display striking differences be- 
tween M 82 and NGC 253, with spectral lines (such as HNCO) 
varying by more than an order of magnitude. We have inferred 
that these differences are related to the heating mechanisms of 
the ISM in the central regions of both galaxies. Our observed 
deficiency of species such as CH3OH, HNCO, NH3, and SiO in 
M82 r elative to other starburst galaxies has been repor ted previ- 
ously dTakano et al.ll2002t IMartm et al.ll2006al l2009al) . In addi- 
tion, we have found a deficiency of NS and HOCO + . The ratio 
NS/CS (proposed as an indicator of shocks in Galactic hot cores, 



IViti et al.ll200"!t) . is found to be an order of magnitude lower in 
M82 than in NGC 253. These abundance differences point to- 
wards a lack of dense molecular material in the central region 
of M82, together with enhanced heating and ionization due to 
UV radiation from newly formed OB stars. These enhanced UV 
fields leave a clear imprint in the molecular composition as il- 
lustrated by the overabundance of species s uch as CO + , HCO , 
HOC + , C 3 H 2 , and CH 3 CCH (reported bv iFuente et alj 120051: 
IMartfn et al]|2006Ul2009btlAladro et al.ll201 lal) It is worth not- 
ing how the abundance of the PDR tracers CO + , HCO, and 
HOC + relative to HCO + (which have been claimed to prove the 
significant presence of PDRs in NGC253, IMartfn et al1l2009bl) 
are similar in both M82 and NGC 253. However, their relative 
abundances in dense molecular gas traced by C 34 S show a sig- 
nificant overabundance in M 82 as a result of the larger amount 
of gas affected by photodissociation. 

The high rotational temperatures measured with NH 2 CN 
transitions in both M82 and NGC 253 are similar to those de- 
rived in Galactic hot cores. If this molecule is tracing the very 
densest gas associated with the hot cores, a higher star-formation 
rate in M 82 relative to that in NGC 253 might be the cause of the 
observed slight enhancement of NH 2 CN in the former. Although 
a significant amount of hot cores are present in M 82, the dense 
molecular gas reservoir flowing into the central few hundred par- 
sees of this galaxy will be exhausted more quickly than in other 
starburst galaxies. However, since NH 2 CN is so faint in M 82, 
further observations of this molecule towards other extragalactic 
sources are needed to confirm its link with dense and hot regions. 

Even if strong differences are found between M82 and 
NGC 253, the comparative s tudy based o n the H NCO/CS ratio 
for a sample of galaxies by IMartm et al.l (l2009al) showed how 
NGC 253 might be in a significantly evolved starburst stage, 
closer to the evolved stage of M 82 than to the younger starbursts 
in galaxies such as M 83. Thus, yet more varied molecular com- 
positions are expected in a wider sample of galaxies including 
young starbursts. 

7. Summary 

We have surveyed the NE molecular lobe of the M 82 star- 
burst galaxy between 129.8 GHz and 175.0GHz and between 
241.0 GHz and 260.0 GHz. In total, we have detected 69 spec- 
tral features corresponding to 18 different molecular species, as 
well as 3 hydrogen recombination lines. The molecules H 2 S, 
CH 3 CN, H 2 CS, NO, and NH 2 CN are detected for the first time 
in this galaxy (for the last two we have tentative detections). In 
addition, we have confirmed previous tentative detections of SO, 
S0 2 , C 34 SandH 13 CN. 

Under the local thermodynamic equilibrium approximation, 
we have derived rotational temperatures and column densities 
of the various molecular species. Although a large number of 
molecules arise from dense and moderately warm gas (n > 
10 4 cirT 3 , r rot = [10 - 30] K), some of them seem to trace the 
cold gas contained in the outer envelopes of the molecular clouds 
(e.g. NO), while others might trace the very warm molecular gas 
phase (e.g. NH 2 CN). 

Making use of the frequency survey of Martin et al. (2006b) 
carried out towards the galaxy NGC 253, we have performed the 
most detailed comparison so far of the chemistry of two star- 
burst galaxies, including twelve species that remain undetected 
in M 82, but for which good upper limits could be obtained. We 
found that, although both galaxies have some similar characteris- 
tics that are typical of starburst galaxies, there are many species 
that point to different dominant processes. For example, C 2 H, 



Aladro et al. (201 1): A A = 1.3 mm and 2 mm molecular line survey towards M 82 



13 



CH 3 CCH, c-C 3 H 2 , or CO + are more abundant in M 82, while 
CH3OH, HNCO, CH 3 NH, NS, HOCO + , and SiO are more abun- 
dant in NGC253. 

We propose that the ratios NS/C 34 S and HOCO + /C 34 S are 
new diagnostics for investigating the evolution of starbursts. 
These molecules, which are not detected in M 82 but are detected 
in NGC 253, could provide valuable constrains on any ultraviolet 
fields or shocks in the ISM of starburst galaxies. 
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Appendix A: Comments on individual molecules 

Eighteen molecular species were detected in our frequency sur- 
vey between 129.8 GHz and 175.0 GHz in the 2 mm atmo- 
spheric window and between 241.0GHz and 260.0GHz in the 
1 .3 mm window. In this appendix, we provide details of the 
Gaussian profile fitting and blending cases of individual tran- 
sitions. 

For some special cases, we performed synthetic Gaussian 
fits using the MASSA software, which reads the spectra and the 
spectroscopic parameters of the selected lines (using the JPL cat- 
alog), and then attemps to measure iteratively the column den- 
sity, excitation temperature, velocity, and FWHM of the lines un- 
til it finds a simultaneous, optimal, Gaussian fits for all of them. 

For instance, the synthetic Gaussian fit method was applied 
to the weak molecules NO, NH 2 CN, SO, and S0 2 . For each 
species, all the transitions that lie in our survey were included, 
enabling MASSA to compare the strengths of the lines and en- 
sure that the line identifications are reliable (since all the transi- 
tions of the same strength are either identified or not). 

In a similar way, synthetic Gaussian fits were applied to 
molecules with hyperfine structure. In these cases, MASSA fits 
one Gaussian profile to each component, taking into account the 
strength of each one. However, since the M82 lines are quite 
broad (60 - 120km s~'), the hyperfine structure is not re solved 
and results in one single Gaussian (see Martin et al. 2010 for de- 
tails about the C 2 H(2 - 1) synthetic fit). This method was also 
applied to these cases where two or more species are blended 
(see Figure lATl for an example). 

The synthetic Gaussian fits do not provide errors in the inte- 
grated intensities. However, these errors are required when do- 
ing the Boltzmann diagrams if one wishes to obtain appropriate 
errors associated with the rotational temperatures and column 
densities. Thus, unless otherwise indicated, in cases where a syn- 
thetic Gaussian profile was fitted, we calculated an error in the 
integrated area using 

(3 x rms x FWHM) / {FWHM/Dv) l/2 , (A.l) 

where rms is the 1 cr noise level of the spectrum that contains the 
line, and Dv is the spectral resolution in velocity units. All the 
parameters are in international system units. 

- Hydrogen recombination lines 

The H36q', H35a, and H34a were observed at 135.3 GHz, 
147.0 GHz, and 160.2 GHz respectively. However, H35a 
is blended with CS (3 - 2) and H36a is blended with 
H 2 CS (4 I>4 - 3i, 3 ). We first fitted a Gaussian profile to H34q-, 
which is the only recombination line that does not suf- 
fer from contamination by other species. We then assumed 
that the three lines have similar intensities, and used the 
same Gaussian profile to subtract H35a and H36a from the 
blended spectra. 

- Cyanoacetylene - HC3N 

We detected five transitions of this linear molecule. The 
HC 3 N(16 - 15) transiti on at 145.561GHz is blended with 
H 2 CO(2(), 2 - lo,i) (Fig. lA.lt . For this cyanoacetylene line, 
we fitted a synthetic Gaussian profile according to the phys- 
ical parameters derived from all the other detected transi- 
tions, and fixed the velocity and the line width to 300 km s 
and 100 km s , respectively. We also included the transition 
HC 3 N ( 24 - 23) at 218.325 GHz observed bv lAladro et all 
(1201 lah . 

- Propyne (methyl acetylene) - CH3CCH 

None of the five detected lines of this symmetric top 



molecule are blended. Each (J — > J -I) transition consists of 
a number of K components (being K = 0, . . . , J — 1) that are 
not resolved because of the large line widths. Thus, we fitted 
unique Gaussian profile to each group of transitions. Given 
the rotational temperature derived from the Boltzmann dia- 
gram (r rot = 31.6 + 0.1K), only the K = to 4 compo- 
nents contribute to the line intensity, with the K = 4 compo- 
nent having a contribution of less than 1% to the total r pea k. 
Therefore, higher values of the /f-ladder (i.e., K > 4) were 
not taken into account. 

Before fitting the Gaussian, we subtracted a bump in the 
CH3CCH(8 - 7) line due to emission from the centre of 
M 82, that was detected by the beam of the telescope. The 
CH 3 CCH(9 - 8) transition, at 153.817 GHz could be con- 
taminated by the HNCO(7oj - 60,6) line, but because no 
other lines of isocyanic acid with similar spectroscopic pa- 
rameters are detected, we assume t hat its contribution is neg- 
ligible. From lAfadro et al.l (1201 lal) we gathered data for the 
CH3CCH (16- 15) transition at 273.420GHz, which lies out- 
side our survey coverage. 

- Carbon monosulfide - CS 

Two transitions of CS fall within our frequency survey. The 
line CS (3-2) at 146.969 GHz is blended with H35o-. We first 
subtracted the estimated contribution of the hydrogen recom- 
bination line, and then fitted a Gaussian profile to the residual 
feature. The CS (5 - 4) at 244.936 GHz does not show any 
special characteristics. To estimate the physical parameters 
of CS, we a lso used the detectio ns reported by iBavet et all 
(I2009bl) andlAladro et al] j20TTah . This is the only molecule 
in the survey that haves two components in the Boltzmann 
diagram, wit h rotational temper atures of 5.8 + 0.1K and 
15.1 + 0.9 K (lAladro et al.ll20lTah . 

- Methyl cyanide - CH 3 CN 

Three transitions of methyl cyanide were clearly detected 
at 147.174GHz, 165.565GHz, and 257.483GHz, none of 
which are blended. In a similar way to propyne, CH3CN is a 
symmetric top molecule whose lines contain a number of un- 
resolved K components (K = 0, 1). As a first approx- 
imation, the Boltzmann diagram gave us T rot ~ 30 K. For this 
temperature only the K = 0, . . . , 4 components contribute to 
the total line intensities, the K > 4 components contributing 
less than 1%. Taking into account only these first /f-ladder 
components, we finally obtain a T mt = 32.9 + 1.6K. 
On the other hand, as mentioned in Sect. [5] this molecule 
probably arises from the cores of the molecular clouds. This 
is consistent with the narrow FWHM of its lines, about 
SOkms- 1 for CH 3 CN(8 JC - 1 K ) and CH 3 CN(9^ - 8 K ), and 
especially for the -60 kms" 1 of the CH 3 CN (14^- 13*:) tran- 
sition. 

- Hydrogen cyanide - H 13 CN 

The group of transitions of H 13 CN(2# - 1#) were detected 
at 172.678 GHz. We fitted a synthetic Gaussian profile using 
the same line width and position for the five hyperfine com- 
ponents of the line. We then assumed r, ot = 20 + 10 K in 
the Boltzmann diagram. The large error in the column den- 
sity reflects the error in the integrated areas found using the 
formula in Eq. Al. 

Following the lower lim it to the carbon iso topic ratio in 
M82, I2 C/ I3 C > 138 bv lMartm etail(l20Tol) . we can esti- 
mate the column density of HC N as A^hcn > 2.3 x 10 14 cm " 2 , 
in agreement with the results of Seaquist & Fraverl J2000). 

- Formaldehyde - H 2 CO 

We detected three transitions of this molecule, the 
H 2 CO(2 u - l u ) line at 140.840GHz being the only that 
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Fig.A.l. Example of a synthetic Gaussian fit performed with 
MASSA. In this case, HC 3 N(16 - 15) (the line on the left) 
is blended with H2CO(2o,2 - 1q,i) (the line to the right). The 
bump at higher frequencies is the emission from the centre of 
the galaxy picked up by the telescope beam. This bump was 
originally subtracted before any Gaussian fit, although here we 
show it as an example of the strength of this emission. See details 
about the Gaussian fitting in the text. 



is not blended. Thus, the width and position of the Gaussian 
fit of this line were used to fix the Gaussian fitting to the other 
two lines. The transition H 2 CO (2 ,2- ln,i) at 145.603 GHz is 
blended with HC3N (16-15), and is also contaminated by the 
emission from the centre of the galaxy detected by the tele- 
scope beam (Fig. IA. lb . We first subtracted the contributions 
of both cyanoacetylene and the emission from the centre of 
the galaxy, and then fitted a Gaussian profile to the residuals. 
Likewise, H 2 CO (2 U - 1 u ) at 150.498 GHz is blended with 
C-C3H2 (22,0 - 1 1,1) and NO (3/2 - l/2)n _ . In this case, we 
first fitted a Gaussian profile to the formaldehyde feature. All 
the fits are synthetic. 

- Cyclopropenylidene - C-C3H2 

A total of eight transitions of this molecule were detected 
in the survey. Six of them are blended with other lines, 
and two are tentative. The transitions c-C3H2(3i,2 - 22,2) 
at 145.090GHz and c-C 3 H 2 (22,o - I1.1) at 150.436GHz are 
blended with CH 3 OH(3 - 2)A+ and H 2 CO(2 u - l u ), re- 
spectively. Both Gaussian profile fits are synthetic. The pair 
of lines at 150.8 GHz is spectrally unresolved. The same 
holds for the pair at 151.3 GHz. In all these cases, the width 
and the velocity of the cyclopropenylidene Gaussian pro- 
files were fixed to 100.0 kms -1 and 300 kms~', respec- 
tively. The fixed velocity is supported by the results of the 
two non-blended lines, at 305.1 and 307.6 kms . However, 
only the C-C3H2 (32,i - 2^2) line at 244.222GHz has a non- 
fixed line width, which is a bit narrower than 100 km s" 1 
(~85kms -1 ). Some molecules occasionally ha ve narrower 
feature s at higher frequency transitions (e.g. CS. lBavet et al.l 
2009b), hence the assumption of a line width of 100km s _1 
may be a good guess for the lower transitions. The tentative 
detections are the two at 151.3 GHz, with line intensities of 
~1 and ~2mK. 

- Methanol - CH3OH 

Six lines or groups of lines of methanol were identi- 
fied at 145.103, 157.179, 157.271, 165.050, 170.060, and 
241.791GHz. A synthetic Gaussian profile was fitted to 
all of them since, except for the lines at 170.060GHz and 
241.791 GHz, all the transitions are tentative, having an in- 
tensity lower than two times the noise level of the corre- 
sponding spectrum. In addition, the CH3OH (3 - 2)A+ line 



at 145.103GHz is blended with c-C 3 H 2 (3i, 2 - 2 2>2 ), whose 
contribution had been previously subtracted. The group of 
CH3OH (Juk-i -Jo,k]E lines at 165.050GHz is also blended 
with the faint feature of SO2 (52,4 -5 1,5). In this case, we first 
estimated the contribution of the methanol lines using their 
spectroscopic parameters and helping us with the other non- 
blended transitions. 

From the Boltzmann diagram, we obtain jVch 3 oh = (1.2 ± 
1 . 1 ) x 1 14 cirT 2 and T^t = 26.2± 1 1 .7 K. The large errors are 
related to the low intensity of the lines because the baselines 
of the spectra have a strong influence on the results of the 
rotational parameters. 

- Nitric oxide - NO 

Four transitions of nitric oxide were tentatively detected 
at 150.176, 150,546, 250.437, and 250.796 GHz. Although 
they are clearly seen (in particular the two at 2 mm wave- 
lengths), their intensities do not reach a S/N of 2. The line 
at 150.546GHz is blended with H 2 CO(2 u - l u ), whose 
contribution was firstly subtracted as explained for formalde- 
hyde. We fitted synthetic Gaussian profiles to the 2 mm and 
1.3 mm hyperfine lines separately. Since we had two groups 
of lines at different frequencies, we were able to calculate 
the rotational temperature. However, this value is not very 
precise because the dynamic range in energies is quite low 
(AE = 5 cm -1 ). Thus, we associated a 50% error to the de- 
rived rotational temperature. The low T Iot ~ 9 K and the high 
column density (N ~ 10 15 cirT 2 ) indicate that nitric oxide 
arises from the external cold layers of the molecular clouds. 

- Thioformaldehyde - H 2 CS 

Only one transition of thioformaldehyde was detected at 
135.297GHz, which is blended with H36<x We first sub- 
tracted the contribution of the recombination line and then 
fitted a Gaussian profile to the residuals. Since we only de- 
tected one transition of this species in the whole survey, we 
fixed the rotational temperature to 20 + 10 K to calculate its 
column density. 

- Ethynyl radical - C 2 H 

Only the C2H (2 - 1 ) transition falls within our frequency 
coverage. It is the brightest spectral feature and the most 
abundant molecule of the survey. This line also shows the 
contribution of emission from the centre of M 82, which was 
firstly subtracted. After that, we fitted a synthetic Gaussian 
profile to the unresolved hyperfine structure, taking into ac- 
count the spectroscopic parameters and line intensities of 
each component. Since we had only one transition of this 
molecule, we assumed T Iot — 20 + 10 K. A stud y dedicated 
to this line in M 82 and NGC 253 is presented in Marti n et ail 
(I20101) . 

- Sulphur monoxide - SO 

Only the SO (43 - 32) transition was detected in the survey at 
138.179GHz. Fixing T mt to 20+ 10K, we obtained a column 
density of 7Y S o = (3.7 + 0.5) x 10 13 cirT 2 . 

- Sulphur dioxide - SO2 

Only the S0 2 (5 2 ,4 - 5 1,5) transition at 165.145 GHz was de- 
tected in the survey. It is blended with a methanol group of 
lines whose contribution was firstly subtracted as explained 
before. After fitting a synthetic Gaussian profile to the resid- 
uals with a fixed position of 300 km s -1 , the SO2 line inten- 
sity was estimated to be ~1 mK. This value is below the noise 
level at this frequency so the detection is considered as ten- 
tative. Because this line is so faint, the Gaussian parameters 
are strongly affected by the baseline. As mentioned at the be- 
ginning of this section, the synthetic fit does not provide any 
error in the integrated area of the line. In this case, we did not 
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use equation Al to calculate the error in the integrated area, 
since the error is even higher than J 7mb dv. Thus, to calcu- 
late the column density we assumed a rotational temperature 
of 20 + 10K, and assigned an error of 50% to the resulting 

- Oxomethyl - HCO 

We detected the HCO (2 , 2 - lo,i) transition at 173.377GHz, 
which contains seven hyperfine lines with lower level ener- 
gies of 2.9 K. It is blended with H 13 CO + (2 - 1). We first 
fitted a synthetic Gaussian profile to the hyperfine lines at 
173.377 and 173.406 GHz because they are less blended with 
H 13 CO + . We fixed the position and line width to 300km s _1 
and lOOkms -1 respectively. In this way, we made sure that 
the synthetic fit covered oxomethyl and not H 13 CO + . We 
then fitted the rest of the hyperfine lines with the same 
Gaussian parameters. To calculate the column density, we 
used a rotational temperature of 20 + 10K. As in the case 
of SO2, the error in the integrated area from equation Al is 
larger than the area itself, thus we did not use any error in 
the Boltzmann diagram. We note, however, that the resulting 
error in the column density is quite large. 

- Oxomethylium, formyl cation - H 13 CO + 

The H 13 CO + (2 - 1) line at 173.507GHz is blended with 
oxomethyl. After fitting a Gaussian profile to the HCO fea- 
ture, we over-fitted another to the residuals, fixing the posi- 
tion and line width to 300kms~' and l OOkms" 1 . U s ing th e 
carbon isotopic ratio in M 82 given by Mar tin et al.l (1201 Ol) . 
12 C/ 13 C > 138, we estimated a column density of HCO + of 
Nnco+ > 1 -8 x 10 14 cm~ 2 . Th i s val ue is consistent with that 
derived bv lSeaquist & Frayed d2000t) . 

- Cyanamide - NH 2 CN 

Five spectral features of cyanamide were tentatively detected 
between 158.815 and 158.943 GHz. The line widths and po- 
sitions in MASSA were fixed to 100.0 and 302.1 km s~\ re- 
spectively, which are the values previously obtained with 
an approximated non-synthetic Gaussian fit. The line in- 
tensities, as calculated with the synthetic Gaussian fitting, 
are about 2 mK. The resulting Boltzmann diagram gave us 
T m = 80.3 + 52.9 K and A? nh ,cn = 1.2±1.5xl0 13 cirT 2 .The 
errors in both parameters, in particular in the column density, 
reflect the strong effect of the baseline on this weakly emit- 
ting species. 

- Hydrogen sulfide - H 2 S 

The H 2 S (1 i,o - 1 o,i ) transition was detected at 168.763 GHz. 
Since there is only one transition of this species, we used 
a r rot = 20 + 10 K in the rotational diagram and derived a 
column density of A^s = (6.1 ± 4.5) x 10 I3 cirT 2 . 
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Fig. A.2. Boltzmann diagrams of the molecules observed in the survey for which we detected more than one transition. The resulting 
rotational temperatures are indicated. 
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0.001 


s 


CH 3 OH(7o^-/-i,x)E 


157.271 


0.12±... 


326.6 ±0.0 


63.2 ±0.0 


0.002 


m, s 


NH 2 CN(8i,7-7o, 7 ) 


159.815 


0.26 ±... 


302.1 ±0.0 


100.0 ±0.0 


0.002 


s 


NH 2 CN(8 1 , 8 -7 U ) 


158.891 


0.26±... 


302.1 ±0.0 


100.0 ±0.0 


0.002 


s 


NH 2 CN(8i,7-7 , 7 ) 


159.943 


0.24 ±... 


302.1 ±0.0 


100.0 ±0.0 


0.002 


m, s 


H34a 


160.211 


1.07 ±0.12 


325.0 ±5.8 


80.0 ±0.0 


0.012 




HC 3 N(18 - 17) 


163.753 


0.80±0.13 


308.5 ±7.5 


87.5 ±14.5 


0.009 




CH 3 OH(./i,^i-./o,k)E 


165.050 


0.13±.. 


330.0 ±0.0 


80.0 ±0.0 


0.002 


m, b, s 


S0 2 (5 2 , 4 -5 u ) 


165.145 


0.09 ±... 


300.0 ±0.0 


78.3 ±0.0 


0.001 


b, s 


CH 3 CN(9 x -8^) 


165.569 


0.48 ±0.14 


268.4 ±12.2 


85.7±31.2 


0.005 


m 


H 2 S (l^o - lo,i) 


168.763 


2.67 ±0.18 


310.4 ±2.7 


80.6 ±6.6 


0.031 




CH 3 OH(3 2 ,i-2 u ) 


170.060 


0.48 ±... 


260.1 ±0.0 


61.3 ±0.0 


0.007 


s 


CH 3 CCH(10^-9^) 


170.906 


9.06 ±0.18 


324.9 ±1.0 


100.0 ±2.4 


0.085 


m 


H 13 CN(2^-U) 


172.678 


0.67 ±0.15 


319.8 ±10.4 


90.3±21.3 


0.007 


hf, s 


HC 3 N(19-18) 


172.849 


0.50 ±0.18 


324.1 ±18.8 


95.8±41.0 


0.005 




HCO(2o, 2 -lo,i) 


173.377 


0.86 ±.. . 


300.0 ±0.0 


100.0 ±0.0 


0.008 


b, hf, s 


H 13 CO + (2- 1) 


173.507 


0.90 ±... 


300.0 ±0.0 


100.0 ±0.0 


0.008 


b, s 


C 2 H(2- 1) 


174.663 


12.11 ± ... 


307.5 ±0.0 


95.6 ±0.0 


0.119 


hf, s 


CH 3 OH(5 ,5-4o, 4 )A+ 


241.791 


0.64 ±... 


343.7 ±0.0 


93.4 ±0.0 


0.006 


m, s 


c-C 3 H 2 (3 2 ,i -2 U ) 


244.222 


0.96 ±0.16 


307.6 ±7.3 


84.5 ±12.6 


0.011 




CS (5 - 4) 


244.936 


4.79 ±0.19 


316.1 ±1.6 


77.1 ±3.6 


0.058 




NO (5/2-3/2)11+ 


250.437 


0.34 ±... 


325.4 ±0.0 


85.1 ±0.0 


0.004 


hf, s 


NO (5/2- 3/2)11- 


250.796 


0.34 ±... 


325.4 ±0.0 


85.1 ±0.0 


0.004 


hf, s 


CH 3 CCH(15k- 14^) 


256.317 


7.18 + 0.13 


312.5 ±0.8 


94.6 ±2.1 


0.071 


m 


CH 3 CN(14*-13*) 


257.483 


0.35 ±0.09 


297.3 ±7.8 


60.3 ±18.0 


0.005 


m 


H 2 CO(2i, 2 -l u ) 


140.840 


0.12 ±0.09 


300.0 ± 0.0 


100.0 ± 0.0 


0.001 


i 


HCO + (3 - 2) 


267.557 


4.18 ±0.27 


278.7 ± 3.2 


102.4 ±7.8 


0.038 


i 



Table A.l. Gaussian parameters for the fit performed to the observed spectral 
features. In those cases where the associated errors are zero, the parameter was 
fixed. Synthetic fits do not give errors in the integrated area. Remarks: (b) blended 
line; (m) multi-transition line. The parameters refer to the main component of the 
group; (s) synthetic Gaussian fit using MASSA; (hf) hyperfine transition; (i) line 
coming from the image band. See Appendix A for more details. 



